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In this short topical review, we revisit a number of works on the pattern-forming dy- 
namical instabilities of Bose-Einstein condensates in one- and two-dimensional settings. 
In particular, we illustrate the trapping conditions that allow the reduction of the three- 
dimensional, mean field description of the condensates (through the Gross-Pitaevskii 
equation) to such lower dimensional settings, as well as to lattice settings. We then go 
on to study the modulational instability in one dimension and the snaking/transverse 
instability in two dimensions as typical examples of long- wavelength perturbations that 
can destabilize the condensates and lead to the formation of patterns of coherent struc- 
tures in them. Trains of solitons in one-dimension and vortex arrays in two-dimensions 
are prototypical examples of the resulting nonlinear waveforms, upon which we briefly 
touch at the end of this review. 



X 



1. Abbreviations 

• BEC: Bose-Einstein Condensate 

• FRM: Feshbach Resonance Management 

• GP: Gross-Pitaevskii (Equation) 

• MI: Modulational Instability 

• NLS: Nonlinear Schrodinger (Equation) 

• OL: Optical lattice 

• RDS: Ring Dark Soliton 

• SI: Snaking Instability 

• TF: Thomas-Fermi 

• TI: Transverse Instability 
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2. Overview 



Bose-Einstcin condensation (BEC) was initially predicted by Bose and Einstein in 
1924 i'2'3 n^ involves the macroscopic occupation of the ground state of a gas of 
bosons, below a critical transition temperature Tc, i.e., a quantum phase transition. 
It took, however, over 70 years for this prediction to be experimentally verified by 
a fascinating series of experiments in 1995 in vapors of rubidium and sodium ^. 
In the same year, first signatures of the occurrence of BEC were also reported in 
vapors of lithium ^ (and were later more systematically confirmed). The ability of 
the experimental groups to controUably cool alkali atoms (currently over 35 groups 
around the world can routinely produce BECs) at sufficiently low temperatures and 
confine them via a combination of magnetic and optical techniques (for a review see 
e.g., ^'^), has been instrumental in this major feat whose significance has already 
been recognized through the 2001 Nobel prize in Physics. 

This development is of particular interest also from a theoretical/mathematical 
standpoint. On the one hand, there is a detailed and "clean" experimental control 
over the produced Bose condensates. On the other hand, perhaps more importantly, 
for experimentally relevant conditions, there exists a very good model partial dif- 
ferential equation (PDE) that can describe, at the mean field level, the behavior of 
the condensate. This mean- field model (which, at heart, approximates a quantum 
many-body interaction with a classical, but nonlinear self-interaction ^'^°) is the 
well-known Gross-Pitaevskii (GP) equation H'l^ rpj^^ latter is a variant of the fa- 
mous Nonlinear Schrodingcr equation (NLS) ^■^, incorporating an external trapping 
potential, and reads: 

z;i*t = -|-A* + 5|'J'|2* + Kxt(r)*, (1) 

where VP is the condensate wavefunction (the atom density is proportional to 
|\['(x, i)P), A is the Laplacian, m is the atomic mass, and the prefactor g, which is 
proportional to the atomic scattering length ^, may take either positive (e.g., for ru- 
bidium and sodium) or negative (e.g., for lithium) values. The two cases correspond 
respectively to repulsive and attractive interactions between the atoms, or to defo- 
cusing and focusing (respectively) Kerr nonlinearities in the language of nonlinear 
optics ^^ . Notice, however, that experimental "wizardry" can even manipulate the 
scattering length using the so-called Feshbach resonances ^^ to achieve any positive 
or negative value of the scattering length (i.e., nonlincarity strength in Eq. (1)). 
Moreover, the external potential can assume different forms. For the "standard" 
magnetic trap, this potential has a typical harmonic form: 

TA '™/22i 22, 2 2n /n\ 

Ve^t = ^{uJ,^x +Uyy +uj^z), (2) 

where, in general, the trap frequencies uJxTUjy,LiJz along the three directions arc dif- 
ferent. As a result, in recent experiments, the shape of the trap (and, hence, the form 
of the condensate itself) can range from isotropic forms to so-called cigar shaped 
traps ^, to quasi two-dimensional (2D) ^^, or even quasi one-dimensional (ID) ^^ 
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forms. Moreover, linear ramps of (gravitational) potential T^xt — rngz have also 
been experimentally applied ^^. Another prominent example of an experimentally 
feasible (and controllable) potential is imposed by a pair of laser beams forming 
a standing wave which generates a periodic potential, the so-called optical lattice 

(OL) 18,19,20,21^ of tjjgfo^j^. 



Vol = Vo 



COs2(-— + 4)^) + COs2(— ii + (t)y) + C0S2(-— 
^x \ ^z 



(3) 



where \x,y.z = Asin(0a;,j/.z/2)/2, A is the laser wavelength, 9 is the (potentially 
variable) angle between the laser beams ^^ and (f>x.y,z are arbitrary phases. Such 
potentials have been realized in one ^^'^^^ two (the so-called egg-carton potential) 
^^ and three dimensions 1^,24 

Moreover, present experimental realizations render feasible/controllable the adi- 
abatic or abrupt displacement of the magnetic or optical lattice trap 22,25 (jjifj^cing 
motion of the condensates), the phase engineering of the condensates in order to 
create vortices ^6,27 qj. dark matter-wave solitons in them 28,29,30,31 ^^^^ stirring (or 
rotation) of the condensates providing angular momentum and creating vortices 
32,33 ^j^^ vortex-lattices 34,35,36^ pj. ^j^g Qiig^YigQ of scattering length (from repul- 
sive to attractive via Feshbach resonances) to produce bright matter-wave solitons 
and soliton trains 37,38,39 goth matter- wave solitons and vortices (which can been 
thought of as the 2D analog of dark solitons) are fundamental nonlinear excitations 
of condensates, and as such have attracted much attention, as can be inferred from 
the relevant experimental and theoretical work. In addition, apart from being inher- 
ently interesting in the BEC physics, these nonlinear excitations are also relevant to 
other fundamental issues of low-temperature (and soft condensed-matter) physics: 
for example we note the dominant role of vortices in the breakdown of superflow 
in Bose fluids 40,4i,42 j^ jg ^^jg^ interesting to mention in this connection that the 
description of such topologically charged nonlinear waves (see e.g., the review *^) 
and their surprisingly ordered and robust lattices, as well as their role in phenom- 
ena as rich and profound as superconductivity and superfluidity were the theme of 
the most recent Nobel prize in Physics in 2003. 

While it is clear that there is an enormous volume of literature on BECs, as well 
as on the nonlinear matter-waves present in this setting, in this topical review, we 
will have to narrow our scope and limit ourselves to the study of "pattern forming 
dynamical instabilities" in BECs, i.e., those that give rise to the generation of 
coherent structures, such as matter-wave solitons and vortex arrays. As we aim 
to focus on the origin of these patterns-which have been observed/realized and 
studied in many different physical contexts as well- we will give a (undoubtedly, 
somewhat subjective) view of the excitement that this new world of "matter waves" 
has brought to diverse disciplines, including atomic physics, nonlinear science, wave 
physics and nonlinear optics, while providing a perspective of the multitude of 
intriguing connections and parallels between these. 

To be more specific, we will first examine the modulation instability (MI) of 
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plane waves, the activation of which, is perhaps the most standard mechanism 
through which bright soUtons and soHtary wave structures appear in the context of 
the "traditional" NLS equation (without the external potential). In this case, the 
plane wave solution of the NLS equation becomes unstable towards the formation 
of a pattern consisting of a chain (or "train" ) of bright solitons. The demonstrations 
of MI span a diverse set of disciplines, ranging from fluid dynamics ^^ (where it is 
usually referred to as the Benjamin- Feir instability) and nonlinear optics *^ (see also 
■^^ for a review) to plasma physics "**. In the BEC context under consideration, we 
aim to present how quasi- ID BECs, trapped in regular magnetic traps or optical 
lattice potentials (in both continuous and discrete settings, as well as in time- 
dependent cases which are relevant to the BEC dynamics in the presence of external 
time-dependent fields), can be destabilized giving rise to such coherent matter- wave 
structures. We will also derive a relevant instability threshold (which will be given 
in terms of physical BEC parameters, thus being quite relevant for experiments), 
providing in this way a simple criterion for avoiding the MI. 

On the other hand, as far as higher-dimensional settings are concerned, the 
most fundamental nonlinear excitations of BECs are vortices, which can be excited 
by phase imprinting, by stirring or rotating the condensate (as mentioned above), 
or by the transverse decay of dark solitons (see, e.g., ^°). The latter mechanism is 
directly connected to the onset of the transverse modulational instability (also called 
simply "transverse instability" (TI) or "snaking instability" (SI)) oi rectilinear darh 
solitons, which has been studied extensively in nonlinear optics (see, e.g., 49,50 ^^^ g^ 
review), as well as in the BEC context ^^. Therefore, the TI is of particular interest 
here, as it is another "pattern forming instability" -relevant in higher-dimension 
BECs-which gives rise to vortex patterns. Here, it is interesting to note that apart 
from the vortex lattices mentioned above, very robust vortex arrays, in the form of 
vortex necklaces composed of vortex-antivortex pairs, can also be formed through 
the onset of TI of ring dark solitons (RDS) in BECs '"'^ (note that relevant patterns, 
but composed by globally linked vortices of one sign rather than of the vortex- 
antivortex type, were also predicted to occur in 2D BECs ^^). In this work, we 
will discuss how the TI destabilizes the RDS, leading to the formation of vortex- 
antivortex necklaces. Also, similarly to the case of the MI in ID settings, we will 
derive the TI threshold (in terms of physical BEC parameters), providing in this 
way a criterion for avoiding the instability. 

The structure of our presentation will be as follows: 

• In section 3, we will discuss reductions of the 3D GP equation to lower dimen- 
sional settings including ID and 2D continuum and lattice ones. 

• In section 4, we will examine the modulational instability in continuum ID set- 
tings and demonstrate how it gives rise to solitary waves. 

• In section 5 we will focus on the modulation instability in discrete settings, and 
in particular in BECs trapped in optical lattices. 

• In section 6, we will turn our attention to 2D scenaria and show how the 
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snaking/transverse instability gives rise to vortex structures. 
• In section 7, after briefly touching upon some of the sohtary wave structures and 
features analyzed in the BEC setting, we summarize our findings and present 
some directions for future study. 

3. Reductions to Lower Dimensions 

The genuinely three-dimensional (3D) BECs can be considered as approximately 
ID if the nonlinear inter-atomic interactions are weak relative to the trapping po- 
tential force in the transverse directions; then, the transverse size of the condensates 
is much smaller than their length, i.e., the BEC is "cigar-shaped" and can be effec- 
tively described by ID models ^'^'^^ Similarly, if the transverse confinement is strong 
along one direction and weak along the others, then for this "pancake-shaped" BEC, 
2D model equations are relevant ^^ . 

Close to zero temperature, it is well-known that the 3D GP equation ^ accurately 
captures the dynamics of the condensate. For the cigar-shaped BEC, the model 
equation is effectively ID and can be expressed in the following dimensionless form: 

In this equation, the normalized macroscopic wave function ij) is con- 
nected to the original order parameter $ through the equation i/; = 
{e/ ^2TT\a\a_\_)~^ exp \iuj±t + i{y'^ + z^)/2a_L] ^, while t and x are respectively mea- 
sured in units of 3/2e^(jjj^ and ^/ia^/2e. Here, ajj_ is the confining frequency in the 
transverse direction, a^ = ^Jh/muji_ is the transverse harmonic-oscillator length 
and e = N\a\/a±_ is a small dimensionless parameter, N being the total num- 
ber of atoms and a the scattering length. Finally, g — ±1 is the renormalized 
nonlinear coefficient, which is positive (negative) for repulsive (attractive) con- 
densates, while V{x) — (1/2)51^2;^ is the magnetic trappping potential, where 
il = (3/2v2e^) • {u)x/^A_) with ojx being the axial confining frequency. Such a 
reduction of the 3D GP equation to ID can be done self-consistently using a multi- 
scale expansion based perturbation technique ^^^^'^ [with the wave function ij) being 
of order 0(1)]; for a rigorous derivation, the interested reader can refer to ^® and 
references therein. 

In the same way, the 2D model for the pancake-shaped condensate assumes the 
form: 

i^^-/^i! + q\i>\'^u + V{r)il} (5) 

ot 

where the roles of the axial frequency lo^ and the spatial variable x are now played 

by the radial frequency wx and r = yjx"^ -\- y"^ respectively; also, the normalized 

variables are connected to the dimensional ones similarly to the ID case, but with 

^x i'^i-) instead of lli±_ {^z)- We note that the use of self-consistent reductions of 

the 3D model to lower-dimensional ones, allows (for different normalizations of the 

wave function) for different choices of the values of coefficients of Eqs. (4) or (5) 
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(e.g., the coefficient of the kinetic energy term may be set equal to f/2, and/or the 
nonUnearity coefficient g may take values different than ±f , and so on; see also a 
relevant discussion in section 6). 

Similar reductions can be performed in the presence of the OL potential: In 
particular, considering an optical lattice along the longitudinal direction of a cigar- 
shaped condensate, the GP model of Eq. (4) is relevant with a potential ^^'^^ 
V{x) = Vb sin (kx), where Vq is the normalized strength of the OL (measured 
in units of the recoil energy Er = /i^/2mA^) and A = 27r/fc is the inter-well distance 
in the interference pattern. In the 2D case, the egg-carton potential of the form 
V{x,y) — Vq [sin (k^x) + sin (fcj,j/)] is the relevant OL potential for the pancake- 
shaped condensates. 

Another useful reduction ^^'^'^ (see also 20,21^ ^^j. periodic (e.g., optical lattice) 
potentials is the one of the GP equation to the discrete NLS equation (DNLS) 
^^'^'^. Starting from Eq. (4) and assuming that the potential V{x) is periodic, i.e., 
V{x + L) — V{x), we consider the associated eigenvalue problem 

, 2° + V{x)(pk^a = Ea{k)ipk.a, (6) 

where (pk.a has Bloch (Floquet) functions (BF's) (pk.a — G'!cp(ikx)ipk.a{x), with 
'4'k,a{x) periodic with period L (a is an index labeling the energy bands Ea(k)). 
For our purposes it is more convenient to use the Wannier functions (WF's) ^^ 
instead of the BF's: 

Wa{x~nL)^J— ipk^a{x)e-'"''^dk. (7) 



271" J-k/L 

Due to completeness of WF's, any solution of Eq. (4) can be expressed in the form 

V'(X, t) =^ Cn,a(t)Wn,aix), (8) 

na 

which after substitution in (4) gives 

' 7, — / ^ '^ni,a^n-ni,a ^ y / , / , ^ni.Qi ^"2,02 ''"3,03 ff 0^10203 ' v^/ 

ni ai ,a2,a3 ni,n2,n3 

where 

/oo 
7/; in in in dr Clfll 

-00 

are overlapping matrix elements. Upon suitable decay of the Fourier coefficients 
and the WF prefactors (which can be systematically checked for given potential 
parameters), the model can then be reduced to 

dc 

* Tr~ — ^O.aCn,a + '-^l.Q {Cn-l,a + Cn+l^a) 

dt 

ai,Q2,a3 
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The latter equation degenerates into the tight-binding model ^°'^^ 
dc 

i ?'" ^ l^O,aCn,a + l^l,a{Cn-l,a + Cn+l,a) + gWJ^222\'^n,a\'^Cn,a, (12) 

if one restricts consideration to the band a only. 

We have now laid out the setting: It can consist of the ID GP equation with a 
magnetic trap or an optical lattice potential (or both); it can also be of the form 
of a 2D GP equation with a magnetic trap and/or an optical lattice. Finally, it can 
also consist of genuinely discrete lattice dynamical systems. In these contexts, we 
now proceed to study some of the pattern forming instabilities. 



4. Modulational Instability in quasi-lD BECs 

4.1. The Case of the Untrapped Condensate 

In the case of the untrapped condensate, i.e., for V{x) = 0, the standard modu- 
lational instability (MI) analysis concerns the investigation of the stability of the 
plane wave solutions of the model equation (4) 45,4647_ j-^ particular, in the absence 
of a potential, Eq. (4) (which is actually the "regular" NLS equation) has solutions 
of the form 'ip{x,t) — ■00 exp[i(fca:; — tot)], with f/^o, k and to being the amplitude, 
wavenumber and frequency of the plane wave respectively, which are connected to 
each other through the dispersion relation lo = k"^ + gtpQ . 

Then, to study the stability of the above mentioned solution, we introduce the 
following ansatz into Eq. (4) 

%p{x, t) = {i]jo + eb) exp[i((fca:; — uut) + ew{x, t))], (13) 

where the amplitude and phase perturbations are assumed to have the form b{x, t) — 
bo ex.p{i{Qx — fit)) and w(x, t) = wq e:Kp{i{Qx — fit)) respectively {Q and fl are the 
perturbation wavenumber and frequency, while bo and wq are constants). Then, 
assuming that the parameter e is small, we can derive [to order 0(e)] the following 
dispersion relation, 

{n-2kQ)'^Q^{Q' + 2gi;^o), (14) 



which signifies that, for g ^ —I, wavenumbers Q < Qcr = "ipoV^^ will be modu- 
lationally unstable, hence the system is unstable to long-wavelength perturbations 
of this form. Importantly, the modulational instability in the absence of potentials 
can only occur in the case of attractive inter-particle interactions {g — —1), while 
it is not feasible for repulsive interactions {g = +1). 

The manifestations of the MI in numerical experiments are shown in Fig. 1. It 
can be clearly observed that MI stable wavenumbers, when perturbed, only sustain 
small oscillations (no sidebands are generated in Fourier space), while unstable 
wavenumbers lead to the generation of soliton trains i.e., arrays of large amplitude 
excitations and the creation of sidebands in momentum space. 
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Fig. 1. The top panel shows the maximum of the solution as a function of time (top left) and 
the Fourier space picture at the end of the numerical experiment (top right) in a modulationally 
unstable case (Q = 1). In contrast, the bottom panels show the same features for a modulationally 
stable case {Q = 2). In both cases, the perturbation 0.05sin((3x) was added to the uniform solution 

V'o = 1- 



Notice that one can take an alternative route to establishing MI, recently intro- 
duced in ®^. Particularly, let us consider the ansatz, 

■0 = [^"0+^(0 cxp(i0a(i)) cxp{iQx)+b{t) exp(i0b(t)) exp(— zQa;)] exp {i{kx — a;i))(15) 

presenting, as previously, a modulation of the plane wave, with a{t), b{t) being 
time-dependent functions. Then, Eq. (15) is introduced not in Eq. (4), but instead 
in the field-theoretic Lagrangian (from which the equation can be derived), namely, 



L 



^(V*Vt-M*)-lV'.P-|9l^l' 



dx. 



(16) 



Deriving the Euler-Lagrange equations and keeping the terms to leading order [i.e., 
to 0{a)], we obtain a = b, (f) = 4>a + <Pb, and 



a = C'2a sin(0) 

(/)=(C1+C3)+2C2COS(0), 



(17) 
(18) 



where Ci = —g'ipl — Q^ — 2Qk, C2 — —gi'o, and C3 = —gipQ — Q^ + 2Qk are constant 
prefactors. Then, one can solve the ensuing ordinary differential equations for a to 
obtain: 



When Q^ + 2gV§ > 0, the solution of Eq. (17) is 




2.9^0 ,• 



y^(Q2 + 2gi^i)QH 



When Q2 + 2^1/^^ < 0, the solution of Eq.(17) is 



«(i) ^ jl - ^ Sinh (^y'(-2.g^2 _ Q2)Q2t^ ' 



(19) 



(20) 
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and hence instability will arise in the latter case (for g — —1 and Q < Qcr = V'oV^), 
while it will be absent in the former one (for Q > Qcr)- It is interesting to note 
that the degenerate critical case of Q = Qcr is also unstable, however the growth 
in the latter is polynomial rather than exponential; in particular, in the latter case 
a{t) ^ ^/l + 4ghp^. We should remark here that this technique can be applied to 
any system exhibiting a Lagrangian structure. 

4.2. The Effect of Linear or Quadratic Potentials 

While the standard analysis of the MI is instructive, in the more realistic case of 
trapped BECs, the interesting variation from the regular NLS equation lies in the 
existence of external potentials in the present context. 

The simplest such example (also one that can be fully analyzed) is the presence 
of a linear (e.g., gravitational ^^) potential V{x) — £x in Eq. (4) ^5,66 jj-^ ^j^jg 
setting, the NLS is well-known to maintain its integrable character ^^. Hence, in 
some sense, we expect that the modulational instability results/conditions will not 
be modified in this case. 

The simplest way to illustrate that is by means of the so-called "Tappert trans- 
formation" ^^ 

ip{x, t) = v{rj, t) exp ( -i£xt - -iS'^t^ ) , (21) 

where 77 = a; + £t^ . This transformation brings Eq. (4) back into the form of the 
regular NLS equation, without the external potential, in which the MI condition 
given in the previous subsection applies. 

The physically more interesting case is the one concerning a quadratic potential, 
say V{x) — /Cx^, with /C being in general a function of time. This case, however, is 
less mathematically tractable. One trick that can be used in this case is the so-called 
lens-type transformation ^^ of the form: 

V'(x,i)=riexp[z/(i)x2]«(C,r), (22) 

where f{t) is a real function of time, C, = x/£{t) and r — T{t). To preserve the 
scaling we choose ^^ 

Tt = l/C". (23) 

To satisfy the resulting equations, we then demand that: 

-/t-4/2+/C(i) (24) 

-IrM- + 4/£2 = 0. (25) 

Taking into account Eq. (23), the last equation can be solved: 

e.{t) = £(0) exp (4 / f{s)ds\ . (26) 
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Upon the above conditions, the equation for v{C,, r) becomes (for g — —1) 

ivr + vqq + \v\^v — 2ijv — 0, (27) 

where f£'^ = 7, and generically 7 is real and depends on time. Thus we retrieve a 
NLS equation with an additional term, which represents either growth (if 7 > 0) 
or dissipation (if 7 < 0) . 

Eq. (27) leads to an explicit, spatially homogeneous solution of the form: 



/ ^0 



if{t)x^ + i 



kx 

W) 



ik^T{t) 



exp 



T{T)+iAli cxp(2r(s))rfs + i6lo 



(28) 



where r(r) = 4 J„ ^{s)ds, and Aq and ^o £^rc arbitrary real constants. 

A particularly interesting case is that of 7 constant. Then, from the system 
of equations (23)-(25) and (27), it follows that /C must have a specific form and, 
thus, the functions /, (, and t can be determined accordingly. In fact, the system 
(23)-(25) and (27) with 7 constant has as its solution JC(t) = (l/l&)(t + t*y^ 
(t* is an arbitrary constant), f(t) = (l/8)(t + i*)"\ i(t) = 2y/2^^/tTt* and 
r(t) — (1/87) In(^i^). This case can be analyzed completely and for details the 
interested reader is referred to 4"^'^^. Importantly, relevant numerical simulations 
corresponding to this case, reveal the formation of solitonic patterns, and in par- 
ticular trains of matter- wave bright solitons, such as the one of Fig. 2 in the case 
of this time decaying potential. 









Ia^ Iaa/ 




Fig. 2. The result of a numerical simulation in a regular (top panel) and semi-log (bottom panel) 
plot for the wavefunction in a potential decaying ^ {t + t*)^'^ . A soliton train is clearly formed. 



In the general case, A must be time dependent (i.e., A = A(t)). Here one cannot 
directly perform the MI analysis and is essentially limited to the realm of numerical 
simulations. Figure 3 shows such an example in the presence of the magnetic trap. 
It can be seen that modulationally unstable wavenumbers result in the formation 
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of soliton wavetrains modulated by the harmonic potcntiaL This is a picture qual- 
itatively reminiscent of the experimental outcome of the modulational instability 
in attractive condensates illustrated in ^"^ . On the other hand, if we initialize the 
system with a modulationally stable configuration, then we note that, interestingly, 
the instability manifests itself in the latter case as well. The reason for the occur- 
rence of MI in both cases is that the dynamics of the potential in Eq. (4) mixes the 
wavenumbers of the original perturbation and eventually results in the excitation 
of modulationally unstable ones. However, this only happens later (because firstly 
the modulationally unstable Q's need to be excited) and with a smaller amplitude 
in this case. 




Fig. 3. The top panel shows the time evolution of the wavcfunction maximum in a (previously) 
modulationally unstable (left) and stable (right) case. The middle panel shows the respective 
spatial profiles at i = 3 and the bottom the corresponding Fourier pictures. Noticeably, the 
instability starts to develop even in the previously modulationally stable case. 



We should note here that the experimental results of '^'^ have significantly in- 
creased the interest in the modulational instability of NLS equations with potentials 
and various additional terms emulating the experimental conditions (see also be- 
low). Such examples can be found e.g., in the works of 69,70,7i,72_ 



4.3. Time-Dependent Settings: Feshbach Resonance Management 

Feshbach resonances 14,73,74 j^g^yg recently become a popular tool for manipulating 
the dynamical evolution of BECs. They have been used to induce or prevent col- 
lapse type phenomena '^3,75^ ^^ study shock waves "^^j and to create condensates 
with tunable interactions (e.g., attractive BECs starting from repulsive ones) more 
generally ^''. The resonance allows one to adjust and essentially tune the scattering 
length [g in the context of Eq. (4)] at will -and even change its sign- by using a 
resonantly tuned ac magnetic field. 

To generate breathing solitons in one spatial dimension (in the presence of the 
magnetic trap 55,78,79^ ^j. jj^ ^j^g quasi-discrete setting created by an optical-lattice 
potential *"), as well to avoid their collapse in two dimensions ^^, an experimen- 
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tally realizable protocol has been recently proposed, in the form of the so-called 
"Feshbach Resonance Management" (FRM). 

The FRM scheme can be modelled (in the mean-field approximation) in the 
framework of the GP equation with the coefficient in front of the nonlinear term 
being a periodic function of time. In Ref. ^*, the periodic function was taken to 
be a piece-wise constant one, periodically jumping between positive and negative 
values. The same model may also be realized in terms of nonlinear optics, where it 
applies to a medium composed of alternating layers with opposite signs of the Kerr 
nonlincarity *^. FRM resembles the dispersion-management (DM) scheme, well- 
known in fiber optics (see, e.g., *^ and ^'^ for relevant reviews), which assumes that 
long fiber links are formed by a concatenation of nonlinear fibers with opposite signs 
of the group- velocity dispersion periodically alternating, thus realizing a system that 
supports robust breathing solitons as information carriers. 

This context motivates the study of the MI in a setting where Eq. (4) has 
time-dependent coefficients. In particular, the most general such setting would be: 

liJt = -D{t)ij,, + 5(i)|^| V + V{x)4,. (29) 

Assuming for simplicity the absence (or weak spatial dependence) of the potential, 
the plane-wave solution to Eq. (29) is 



ipo = ^0 exp 



i[kx- k'^ / D{s)ds ~ Al gis)ds 
Jo Jo 



(30) 



Note that the arbitrary amplitude Aq in Eq. (30) can be rescaled to ^o = 1- Then, 
we seek, as in the previously studied cases, solutions of Eq. (29) incorporating a 
modulation that perturbs the above mentioned plane wave, having now the form 

Vj = -00 [1 + ew{t) cos(Qa;)] , 

where e and Q are the amplitude and wavenumber of the perturbation. This way, 
we derive the following (linearized) equation for the real part, Wr = Re(w), of the 
perturbation, 

Wr = DD-^Wr - Q^D{t) [Q^D{t) + 2g{t)] Wr, (32) 

where the overdot stands for d/dt. There are several special cases of this equation 
that were previously studied. In the case of DM (i.e., for D = D{t) and a{t) = 
const.), the MI analysis was performed in ^^. On the other hand, in the FRM 
context, for D = 1/2 and time-periodic g{t), Eq. (32) is a Hill equation that was 
considered in ^*, while the more specialized case oi g{t) = 1 + 2acos(wi) leads to 
the Mathicu equation that was studied, in this context, in ^^. 

Motivated by the FRM scheme proposed in ^*, Ref. *^ discusses the special case 
of a piecewise-constant time-dependent g alternating between gi (for < t < r) 
and —52 (for r < i < T) as mentioned above. In this case, first, we define sf = 
Q'^iQ'^/4: + gi) and s\ = Q^(g2 - Q^/4), which assumes Q^ < 4g2- Then, seeking 
for a solution in accordance with Bloch's theorem ^'^ within one period and using 
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the appropriate continuity conditions for the integration constants, one can obtain 
a stabihty condition. This procedure (similar to the one used for identifying the 
ahowed/forbidden bands for the Kronig-Penney model in quantum mechanics *^) 
leads to the determinant condition for the eigenfrequency (Floquet multiplier) fi of 
the solution: 



cos(rJT) 



2siS2 



■sin(sir)sinh[s2(T-T)]+cos(sir)cosh[s2(r-r)] = F(Q).(33) 



If the above condition Q^ < 452 does not hold, we redefine S2 — \/Q^{Q'^/4: — 52), 
and obtain, instead of Eq. (33), 



cos{nT) 



S2 



2siS2 



sin(siT) sin[S2(T - r)] + cos(sir) cos[S2(T - r)] = F(Q). (34) 



By examining the function |i^((3)| or 



FiQ) 



defined in Eqs. (33) and (34), and 



comparing it to 1, we can find whether there is a real eigenfrequency for a given 
perturbation wavenumber Q, or it belongs to a "forbidden zone" , which implies the 
MI. 

A typical example of results produced by this analysis is shown in Fig. 4, for 
91 — 92 — 0.3 and r = T/4 = 1. In general, we have found that the number 
and widths of the instability windows increase as long as the mean value g = 
[giT — 52(7 — t)]/T of the scattering length gets large and negative, i.e., when the 
BEC is "very attractive" on the average. 




Fig. 4. The plot of the function F{Q) defined by Eqs. (33) and (34) for gi = 32 = 0.3 and 
r = T/4 = 1. When the function satisfies |i^(Q)| < 1, the perturbation of wavenumber Q is 
modulationally stable, while for [_F(Q)| > 1, it is modulationally unstable. 



In *®, the results of the theoretical predictions for the MI were tested against 
direct numerical simulations of the original PDE and the agreement between the two 
was found to be generic. Furthermore, the case of the presence of a weak magnetic 
trap was also examined. 



February 2, 2008 8:52 mplb6a 



14 Kevrekidis and Frantzeskakis 

4:A. Avoiding the Modulational Instability 

In all of the above settings, we examined the conditions for inducing a modulational 
instability, in order to form solitary wave patterns due to its presence. However, 
in other contexts, it may be desirable to avoid such dynamical manifestations and 
have robust structures that are not subject to MI. In such cases, a simple criterion, 
involving physical parameters of a trapped condensate, for the suppression (or the 
onset) of the MI would be useful, especially for relevant experiments. As proposed 
recently in *^, one can derive such a criterion in terms of the magnetic trap strength 
(which controls the size and the number of atoms of the condensate) . 

In particular, as discussed in section 4.1, the modulational instability sets in for 
perturbation wavenumbers below a critical one Q < Qcr = v2'0o [we set g — —1 
in the framework of Eq. (4)]. This can be equivalently interpreted as follows: when 
length scales 

\>X„ = y^ (35) 

become "available" to the system, then the modulation over these scales leads the 
plane wave solution to instability. 

However, in the presence of the magnetic trap, there is a characteristic scale 
set by the trap, namely the BEC axial size, XbeCj which depends on the trap- 
ping frequency fJ. When Xbec < ^cr, suppression of the modulational instability 
is expected. To estimate Xbec in a specific setup, one can examine e.g., a pro- 
tocol relevant to the recent experiments reported in ^^ (to which we referred ear- 
lier). In particular, the initial configuration consists of a ID repulsive condensate, 
whose ground state wave function can be derived in the framework of the so-called 
Thomas- Fermi (TF) approximation ''; subsequently the interaction is made attrac- 
tive, by taking advantage of the Feshbach resonance. Therefore, in this situation 
the length scale of the condensate can be calculated in the TF approximation as 
Xbec ~ 2^/2/7/17, where /i is the chemical potential. Therefore, if 

n>ncr^ — K 0.387 (36) 

(where we took into regard that /i — (3/2)^' "^il^'"^, resulting from the normalization 
condition of the wave function ip) , then the trapping conditions are "engineered" in 
such a way that the modulational instability cannot manifest itself. Such stringent 
trapping conditions were indeed found to preclude the occurrence of the MI in the 
numerical experiments of **. 

5. Modulational Instability in BECs trapped in Optical Lattices 

5.1. Discrete Modulational Instability Theory and Experiment 

Another interesting manifestation of the modulational instability has recently been 
theoretically predicted and experimentally verified in the context of BECs trapped 
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in an optical lattice. The latter setting, as discussed in the section 3, can be studied 
(under appropriate conditions) in the framework of the tight binding model of Eq. 
(12). Hence, we will examine the instability in the genuinely discrete context to also 
illustrate how the relevant conditions are modified by discreteness. 

The relevant discrete reduction of the equation can be re-written in the general 
form: 

^^^ = "^(V'n-l + V-n+l) + (K + .g I V'« l')V'n, (37) 

at 
where ^„ is the normalized wavefunction of the condensate in the n-th site of 
the array, while the constant D describes the macroscopic tunneling rate between 
adjacent sites. In the context of Eq. (37), for weak harmonic potentials (such that 
we can consider Vn — Q for the fraction of the lattice considered) the MI analysis can 
be carried out ^^ similarly to the continuum case. In particular, taking into regard 
that the plane wave solution of Eq. (37) has the form ■(/;„ — ipQ exp[i(fcn — wi)] (with 
u! = —2D cos{k) + giV'oP) and performing the stability analysis along the lines 
of section 4.1 with a perturbation ~ exp(z(5n), the following condition is readily 
derived: 



{n ~ 2D sm{k) sm{Q)y = 8L»cos(fc) sin2(-^) 



2i?cos(fc)sin2(|)+g|V;oP 



(38) 



The crucial difference of this equation from its continuum analog (cf. Eq. (14)), 
identified in ^^ and in the context of BECs later in ^^ , was the existence of unstable 
wavenumbers, even for the case oi g > i.e., for repulsive condensates. In partic- 
ular, it was appreciated that should wavenumbers tt/2 < k < 3tt/2 be excited in 
the repulsive system, it will then have the potential for instability provided that 
sIV'oP > — 2Z?cos(fc) sin^(g/2) (i.e., for a sufficiently large amplitude excitation). 
Hence the key to the instability is to excite the modulationally unstable wavenum- 
bers larger than tt/2. Since the experiments are conducted in the presence of a 
magnetic trap, an indirect way to excite such wavenumbers was suggested in ^"^ and 
subsequently experimentally implemented in ^^. In particular, if the magnetic trap 
is instantaneously displaced, then the BEC will have to move towards the new cen- 
ter of the magnetic trap. This will excite a quasi- momentum (i.e., a wavenumber) 
which is small if the displacement of the center of the trap is small, but can become 
larger than tt/2, for sufhciently large initial displacements. By using the Josephson 
equations 

h—^{t) =21? sinA<^(t) (39) 

hj^Acj^it) = -2n m, (40) 

for the collective coordinates of ^ = ^ ■ jrij (center of mass) and <i>j+i (t) — (pj (t) = 
A(/)(t) (quasi- momentum) , where fl — (l/2)m(x'^(7r/fc)^ is the relevant trap fre- 
quency, one can analyze the BEC motion. In particular, it can be obtained ^° that 
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^>Ccr = \ -FT- 



(41) 



in order for A0 to becorae 7r/2, a displacement of 

~n 

should be applied to the center of the trap. This theoretical prediction was tested in 
the numerical experiments of ^°, both in the discrete equation and in the continuum 
one with the periodic potential, yielding very good agreement with the numerical 
findings (for details see Fig. 5). The theoretical results of ^° were also tested directly 



0.04 


(a) 




let 










\ \ 


J 


111 


l>r^:\ 


\ 


(b) 




w 




OJ 












M 


m. 


j) 


L.K'^% 



^^-n "'^ 


/ 


Ilk 


\ 


90 /^l 


40 




\ 








so 





-20 I) 111 



111 4i> m 

X (fUll) 



100 

time (ms) 



Fig. 5. Left panel: the density \ipn\'^ at different times 0,20,40 ms (from the right to the left of 
the figure) for initial displacements of (a) 80 and (b) 100 sites, respectively below and above the 
critical displacement; Right panel: plot of the quasi-momentum < A: > vs. time for three different 
initial displacements: 40, 80 and 90 sites. The critical displacement is ^cr ~ 84 sites. It is clearly 
seen that when < fc > reaches 7r/2 (i.e., for an initial displacement greater than ^cr), the MI is 
activated. 



in experiments of the Florence group, published in ^^. The comparison and good 
agreement between the experimental results and the theoretical predictions can be 
summarized in Fig. 6. 

We should also note that a different approach to the modulational instability in 
optical lattices was recently proposed in ^^ , where it was observed that the merging 
of the regular plane wave branch with an unstable, period-doubled branch of states 
is responsible for the onset of the instability. It would be interesting to examine 
whether this type of scenario can carry over to different settings. 



5.2. Modulational Instability in a Time- Dependent Optical Lattice 

The context of the optical lattice lends itself to an interesting variant of the FRM 
scheme discussed above. In particular, for example, the heights of the interwell bar- 
riers of the optical lattice are proportional to the intensity of the lasers and can be 
easily modulated periodically in time. This induces an oscillating tunneling ampli- 
tude of the condensates between adjacent wells, as well as an oscillating interaction 
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Fig. 6. Comparison of experimental results with Eq. (41). Filled circles represent coherent oscil- 
lations, empty circles denote pinned motions; continuous line results illustrate the critical line of 
Eq. (41). 



energy of the condensates trapped in each well. This possibility along with the more 
conventional FRM scheme presented above motivated the recent study of ^^ where 
the coefficients of Eq. (37) were allowed to vary in time according to the form: 

-i:'(t)(V„+l + V-n-l) + gWllZ-nlV™, 



ill>n 



(42) 



Since one of the two time dependences can be absorbed in a rescaling of time, Eq. 
(42) was examined in ^^ with Dit) = 1; additionally, in the same work, a periodic 
variation of the form g(t) = 1 + esin(a;t) was assumed for the time-dependent 
nonlinearity coefficient g{t). 

To examine MI in this context, we first note that the plane wave solution of Eq. 
(42) reads: 

V (t) = g«('="+2 cos(fe)i- Jj,'g(s)ds)^ 



Then, we consider a perturbation ansatz of the following form: 
Mt) = vn{t) [1 + ~e (a(t)e'«" + /3(t)e-^Q")] , 



(43) 



(44) 



where Q is the perturbation wavenumber and a{t), j3{t) are complex, time depen- 
dent fields. The resulting equation for a is then of the form: 

a -f {2iB — euj cos(cji)) a 
+ [A^ -B'^ + 2A + e{2A sin(cjt) ~iuj{A + B) cos(t^i)] a = . (45) 

where A — 2 cos fc(l — cos Q) and B ^ 2 sin k sin Q are constants depending only on 
the wavenumber k and the perturbation wavenumber Q. To examine this periodic 
coefficient equation (for possible parametric resonances), a multiple scales technique 
^^ was implemented in ^^, using k — kg + efci + O(e^). It was thus found that 



-1 



fci = ±- 



i + ic.2 



4,/l -I- 1:^2(1 -cosQ)sinfco 



(46) 
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Therefore, the boundaries for parametric instabihties on the (fc, e) plane can be 
identified as: 



4i/l + 7^2(1 - cos(5)sinfco 
±e - -^ . (fc - fco), 



-f 



f 



\-' 



(47) 



and the resuh is demonstrated in diagrammatic representation of Fig. 7. The points 
in the figure represent examples of wavenumbers in the different regions of parame- 
ter space. In ^^ the numerical results indicated good agreement with the theoretical 
prediction except for the region between the parametrically and modulationally 
unstable intervals (where higher-order parametric resonances may set in). Finally, 
the difference between the manifestations of the MI and of the parametric insta- 
bility which can occur in this setting, is that the former has a considerably larger 
growth rate and hence develops much earlier (as was analytically calculated and 
numerically verified). 




Fig. 7. The plane of (e, k) where the region of the parametric instabihty represented by Eq. (47) is 
shown by two straight hnes. Wavenumbers between the two "stability limits" will be (to leading 
order) parametrically unstable. 



6. Transverse (Snaking) Instability in quasi-2D BECs 

6.1. Transverse Instability of Dark Solitons 

We now turn our attention to higher-dimensional settings and discuss a similar, 
long-wavelength, instability in the 2D setting, relevant to repulsive pancake-shaped 
BECs. In particular, we will consider the transverse modulational instability of dark 
solitons, which results to a transverse snake deformation of the soliton (hence the 
name "snaking instability" (SI)) ^^'^o, causing the nodal plane to decay into vortex 
pairs. 

Before proceeding to demonstrate how the SI of dark solitons may lead to robust 
vortex patterns (vortex necklaces) in BECs, it is relevant to recall some qualitative 
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results concerning the SI. First, it is convenient to renormalize Eq. (5) (for g = +1) 
to put it into the form ^^, 

dib 1 

i^^--l^ij + \M''u + V{r)i> (48) 

In the case of the untrapped condensate, and particularly on a homogeneous back- 
ground density f/'oi ^'^ exact analytical solution of Eq. (48) in the form of a recti- 
linear dark soliton along the x-direction reads ^^, 

V'(a;, t) — "00 (cos(^ tanhC + i snap) cxp{—ifit), (49) 

where fj, = ipQ is the chemical potential, (p is the phase shift {\(p\ < 7r/2) across the 
dark soliton, (^ = ipo (cos(/3) [x — i/'o(sin(^)i], while the amplitude and velocity of the 
soliton arc given by cos ip and sin ip respectively; note that the limit case </? = 
corresponds to a stationary dark soliton, u = V'o tanh(-!/)oa;) exp(— z/it), often called 
"black soliton" . In the same case (i.e., in the absence of the potential), the transverse 
modulational instability of the dark soliton in Eq. (48) occurs for perturbation 
wavenumbers 



< 



>y sin^ (/> + -(/iq 2 sin (/) + Vo "^ - (2 sin (/> + Vq ^) 



-,1/2 



(50) 



Equation (50) shows that the instability band may be suppressed if a dark stripe 
is bent so as to form a ring of length L < 2itQ. This argument led to the concept 
of ring dark solitons (RDS) ^^, which (under certain conditions) are not subject to 
the SI. The RDSs were first introduced in the context of optics, where they have 
been studied both theoretically ^® and experimentally ^^, while, later, they were 
also predicted to occur in BECs ^^. Importantly, in the latter context, the presence 
of the magnetic trapping potential is responsible for the formation of robust vortex 
necklaces, when the RDS are subject to the snaking instability. 

6.2. Vortex Necklaces formed due to the Snaking Instability 

A ring dark soliton solution of Eq. (48) can be sought in the form ip = ipTpv ^^, 
where 



4'TF = Va^^(i72H^^cxp {-^^lt) , (51) 

is the approximate ground state of the system (Thomas-Fermi cloud) and fi is the 
chemical potential; also, consistently with the form of the rectilinear dark soliton 
in Eq. (49)), v is the wave function of the RDS, given by 

v{r,t) = cos ip{t) ■ tanhC + ism(p{t), (52) 

where C = cos(p{t) [r — R{t)], while p{t) and R(t) arc the time-dependent phase 
and radius of the RDS. Then, the adiabatic perturbation theory for dark solitons 



February 2, 2008 8:52 



mplb6a 



20 Kevrekidis and Frantzeskakis 

''^'^^ can be used to obtain (for n = 1 and almost black solitons with cosi,^ « 1), 
the following equation of motion for the time-dependent ring radius, 

d^R _ I dWjR) , 1 

dt^ " 2 dR +3i?- ^ ^ 

This equation implies that the RDS moves in a combined effective potential ^^ 
U{R) = (l/2)(fii?)2 - (1/3) Ini?. Furthermore, there is a critical radius Re = 
■\/2/(3il), for which the RDS in the BEC is stationary (unlike its optical coun- 
terpart which expands indefinitely ^^); this radius can also be obtained through 
solvability conditions ^^. For any R y^ Re, the ring will execute oscillations between 
a minimum and a maximum radius in the potential n(_R) ^^, with the frequency of 
small oscillations around the minimum being Q. 

However, while for shallow rings (initially placed near the BEC center) the 
dynamics follows the above effective description (cf. left panel of Fig. 8), for suf- 
ficiently deep rings, the RDS dynamics is much more complex. In particular, the 
rings become unstable to azimuthal perturbations and, thus, subject to the SI that 
results in their breakup. This, in turn, leads to the formation of vortex-antivortex 
pairs. How many such pairs will be formed originally is a function of the soliton 
depth. In ^^, it was found that 8, 16, 24 or 32 vortices could initially be formed, 
from the RDS breakup, but eventually only 8 vortex patterns would be selected by 
the long-time asymptotics of the system (cf. right panel of Fig. 8). The latter were 
found to dynamically alternate between elaborate x- and cross-shaped patterns ^'^ . 
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Fig. 8. Left panel: Evolution of a moving (shallow) ring dark soliton. The top left subplot shows a 
cross-section of the radial dependence of the wavefunction (TF+RDS). The top right, bottom left 
and bottom right subplots show the evolution of a ring, as it shrinks and subsequently expands. 
Right Panel: Evolution of a deep, originally stationary ring soliton, breaking up due to the snaking 
instability into (eventually) 4 vortex pairs. 
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6.3. Avoiding the Transverse Instability 

Similarly to the case of the modulational instability in quasi-lD condensates (cf. 
section 4.4), in this subsection, we aim to illustrate the criterion for suppressing 
the transverse instability of dark solitons in quasi-2D BECs (see relevant work in 
^^ and from that perspective also in **). Furthermore, we will test this criterion 
against direct numerical simulations, exposing the possible dynamical scenarios and 
quantifying their dependence on the trapping parameters. 

To proceed in this direction, we confine ourselves to the case of stationary (black) 
solitons with sincj) — 0, for which Eq. (50) leads to Qcr = ipQ^ (recall that ipQ is 
the density of an untrappcd BEC). On the other hand, in the case of a harmonic 
trapping potential V{r) = (l/2)0^r^, a similar calculation as for the ID problem 
yields the characteristic length scale of the BEC (i.e., the diameter of the TF cloud) 
as Xbec ~ 2y^2fl/n. Then, the criterion for the suppression of the transverse 
instability is that the scale of the BEC is shorter than the minimal one for the 
instability. The corresponding condition reads 



n>^. (54) 

To obtain the minimum value of Q we need to know how the chemical potential /i 
is connected with ipQ. As a first guess, in the absence of the dark soliton, one can 
assume V'o ~ M (close to the center of the BEC), which yields ft > V2/tt = 0.45. 
Hence, stronger trapping should "drown" the transverse instability and preserve 
dark soliton stripes on top of the Thomas- Fermi cloud (i.e., stable "dipole" solu- 
tions). Note that in terms of real physical units, the above mentioned critical value 
of fl may correspond, e.g., to a weakly interacting ^^Rb pancake condensate, con- 
taining w 10'^ atoms, confined in a trap with w^ = 27r x 5Hz and w^ = 27r x 50Hz. 
This condition was numerically tested in *^ and was found to be an overestimate of 
the critical trapping frequency for the transverse instability, of f2cr ~ 0.31. It was 
also found in ** that for 0.18 < fi < 0.31, while the stripe was dynamically unsta- 
ble, there is not sufficient space for the instability-induced vortices to fully develop; 
as a result, after their formation, they subsequently recombine and disappear (cf. 
bottom panels of Fig. 9). On the contrary, for O < 0.18, the vortices will survive 
in the asymptotic evolution of the system, and, naturally, the weaker the trapping 
the larger the number of "engulfed" vortices generated due to the stripe breakup 
(cf. top right panel of Fig. 9). A question that naturally arises in the results above 
concerns the disparity between the critical point theoretical estimate for the trans- 
verse instability and the corresponding numerical finding. In ^®, this disparity was 
attributed to the fact that the theoretical stability analysis of ^^ is performed for 
the infinite homogeneous medium, in the absence of a magnetic trap. In the pres- 
ence of the trap on the one hand, the background is inhomogeneous, while, on the 
other hand, for tight traps resulting in small condensate sizes, the presence of the 
dark soliton at the BEC center, significantly modifies the maximum density. Thus, 
one should expect that the relation ipQ = fi should be modified as SipQ = M; where 
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the "rescaling" factor 6 < 1. The numerical simulations indicated that embedding 
the dark soliton in the BEC center reduces the maximum density of the TF cloud 
to half of its initial value. This suggests that 5 = 1/2, which, in turn, leads to a 
new minimum value of fi, namely f2 = I/tt « 0.318. This modified criterion for 
the suppression of the transverse instability is in a very good agreement with the 
numerically obtained condition of fJ > 0.31. 




Fig. 9. Top: the panels show the contour plots of the density |i/ip for n = 0.35 at t = 1000 (left) 
and n = 0.15 at i = 200 (right). In the first case the transverse instability is clearly suppressed, 
while in the second it sets in, giving rise to a formation of a vortex pair. Bottom: Snapshots of 
a vortex-pair evolution in a case where snaking instability has set in (H = 0.2). In the left panel 
(t = 190) the formed vortex pair is shown, while the right panel (t = 210) shows the recombination 
of the two vortices, resulting in the re-creation of a dark stripe structure. The latter is unstable 
and decays at longer times (t fs 400). 



7. Outlook 

In this short review, we have tried to touch upon a number of aspects of certain 
( "pattern forming" ) dynamical instabilities that Bose-Einstein condensates can sus- 
tain. The main message that we have tried to convey is that the mean field models 
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of BECs result into nonlinear partial differential equations of the NLS type, about 
which a wealth of features is inherited from earlier studies in nonlinear optics and 
nonlinear wave theory. However, this new exciting world of matter-waves also offers 
a number of interesting twists (such as the default inclusion of trapping poten- 
tials, or the external control via magnetic fields) that can be interwoven and may 
considerably modify the previously developed body of knowledge. 

A particularly exciting feature that we chiefly tried to explore in this novel con- 
text of BECs in the present work, is the potential of dynamical transitions that, in 
turn, lead to the formation of patterns. Such patterns may range from matter- wave 
bright soliton trains and bound states in quasi-one dimensional systems to vortex 
compounds and vortex lattices in higher-dimensional settings. The experimental 
realization of such structures 34,35,36,37,39,91 j-^^^g j-g^Hy motivated and intensified 
the detailed examination of nonlinear waves in these contexts. Matter-wave bright 
^'^ and dark i"0'ioi solitons, twisted sohtons 102,103^ g^^p goljtons i04,i05^ Feshbach 
solitons '^^''^9, propeller domain-wall solitons ^°^, shock waves i07, 108,109^ vortices in 
the presence 110,111 ^j, -j-j^g absence n^^ns of optical lattices, interactions of solitons 
with sound ^^^ or with other solitons n^'H^ are only a few examples of structures 
that have been studied in this context. A useful resource letter containing many 
relevant references (especially for vortices and various excitation modes of BECs) 
can be found in ^^'^ . 

In conclusion, nonlinear waves and patterns have really found a new paradigm 
in the context of Bose-Einstein condensates. Yet, numerous aspects of this general 
direction still need to be clarified. Such examples include (among many others): 

• aspects of the interaction of the condensate with the gas for nonzero temperatures 

118,119,120. 

• features of the collapse type phenomena that can occur in higher dimensional 
settings (where experimental results still challenge the theoretical interpretations 

73). 

• multi-component condensates where numerous novel structures still await to be 
discovered/observed 121,122,123^ 

• as well as fundamental instabilities such as the Landau instability, that have 
sparked a long standing controversy that still waits to be undisputably resolved 
(see e.g., the recent discusion of ^^4 ^^(j references therein). 



This exciting journey has just begun 
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Note Added in Proof 

After the submission of this paper, more detailed experimental results illustrat- 
ing the appearance of the modulational instability in a ID moving optical lattice, 
in good agreement with the Gross-Pitaevskii theory, were presented in ^^^. 
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